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Abstract 

We construct a representation of the coherent state path integral using the Weyl symbol of the 
Hamiltonian operator. This representation is very different from the usual path integral forms 
suggested by Klauder and Skagerstan in 1], which involve the normal or the antinormal ordering 
of the Hamiltonian. These different representations, although equivalent quantum mechanically, 
lead to different semiclassical limits. We show that the semiclassical limit of the coherent state 
propagator in Weyl representation is involves classical trajectories that are independent on the 
coherent states width. This propagator is also free from the phase corrections found in |2| for 
the two Klauder forms and provides an explicit connection between the Wigner and the Husimi 
representations of the evolution operator. 

PACS numbers: 03.65.Db, 03.65.Sq 
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I. INTRODUCTION 



The set of coherent states forms a non-orthogonal over-complete basis. This has im- 
portant consequences for the path integral formulation of the propagator. It implies the 
existence of several forms of path integrals, all equivalent quantum mechanically, but each 
leading to a slightly different semiclassical limit. Klauder and Skagerstam (KS) l| proposed 
two basic forms for the coherent state path integral, whose semiclassical limits were consid- 
ered in jjj. It was shown in Q that these two semiclassical propagators can written in terms 
of classical complex trajectories, each governed by a different classical representation of the 
Hamiltonian operator H: the P representation Hp in one case and the Q representation 
Hq in the other. We briefly review these representations and their semiclassical limits in 
section 2. The two most important characteristics of these semiclassical formulas are, first, 
that the underlying classical dynamics depends explicitly on the width of coherent states. 
Second, the phase appearing in these semiclassical formulas is not just the action of the 
corresponding complex classical trajectory, but it also contains a 'correction term' / that 
comes with different signs in each formula (see Eqs. (jl5J) and 

In (2 1 it was also suggested that a semiclassical representation involving directly the Weyl 
representation of H, or the classical Hamiltonian Hyy, could probably be constructed, and 
a formula for this representation was conjectured. A first attempt to derive such formula 
was recently presented in jjj. The strategy used there was to build the propagator out of 
infinitesimal propagators that alternated between the two KS forms. The resulting semiclas- 
sical dynamics turned out to be governed by (Hq + Hp)/ '2, which coincides with the Weyl 
symbol for polynomial Hamiltonians with up to cubic terms in q and p only. The correction 
to the action was found to be (Iq — Ip)/2, which is also non-zero for general Hamiltonians. 
In this paper we construct a new representation of the quantum mechanical path integral 
in the coherent state representation that contains precisely Hw and derive its semiclassical 
imit. The new construction is based on the properties of translation and reflection operators 

in 

% |5[ , which form basis for expressing general operators. While in the KS path integrals each 
path contributes a term of the form expiS/h, where S is the action along the path (com- 
puted with either Hq or Hp), the exponent in the new form is rather different and does not 
immediately resembles an action. Although the terms in this exponent can be re-arranged 
so as to look similar to the action function, it is only when the limit of continuous paths is 
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sentation is indeed given by the expression conjectured in |2j: the underlying dynamics is 
purely classical (independent on the width of the coherent states) and there is no correction 
term to be added to the action. More importantly, the new path integral representation 
allows for a direct connection between the coherent state representation of the evolution 
operator and its Weyl symbol. 

The paper is organized as follows: in section 2 we review the path integral constructions 
of Klauder and Skagerstan and their semiclassical approximations. In section 3 we construct 
the new path integral representation and in section IV we derive its semiclassical limit. The 
two path integrals of Klauder and Skagerstan are compared with the new form in section 
V, where we also comment on relevance of these results for numerical calculations. Finally, 
in section VI, we discuss the connection between the Weyl symbol of the evolution operator 
and the diagonal coherent state propagator. 

II. THE COHERENT STATE PROPAGATOR AND ITS SEMICLASSICAL AP- 
PROXIMATIONS 

In this section we define the coherent state propagator and review the construction of the 
two path integrals suggested by Klauder and Skagerstan, showing how the symbols Hq and 
Hp of the operator H appear in each of them. We also write down the semiclassical limit of 
these path integrals to compare with our results in the next section. Our presentation here 
is strongly based in (^j]. 

A. The propagator 

The coherent state \z) of a harmonic oscillator of mass m and frequency uj is defined by 



^) = e -il"l 2 e" at |0) 



(1) 



with |0) the harmonic oscillator ground state and 
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In the above q, p, and er are operators; q and p are real numbers; z is complex. The param- 

i i 

eters b = (h/mu;) 2 and c = (hmu) 2 define the length and momentum scales, respectively, 
and their product is h. 

For a time-independent Hamiltonian operator H, the propagator in the coherent states 
representation is the matrix element of the evolution operator between the states \z') and 
\z"): 

K{z",z',T) = {z"\e-^ T \z'). (3) 

We restrict ourselves to Hamiltonians that can be expanded in a power series of the creation 
and annihilator operators <v and a. 

In the construction of a path integral for K, and also in the derivation the semiclassical 
limit of the propagator, the Hamiltonian operator H is somehow replaced by a classical 
Hamiltonian function H(q,p). This 'replacement', however, is not uniquely defined, and the 
ambiguities that exist in the relation between the operator H and the function H(q,p) also 
arise in connection with the overcompleteness of the coherent state basis, as we shall see in 
the next subsections. 

There are actually many ways to associate a classical function of position and momentum 
A(q,p) to a quantum mechanical operator A p. However, three of them are specially 
important. The first one, denoted Ag(q,p) and called the Q representation of the operator 
A, is constructed as follows: one writes A in terms of the creation and annihilation operators 
a' and a in such a way that all the creation operators appear to the left of the annihilation 
operators, making each monomial of A look like c nm a) n a m . Then we replace a by z and 
a) by z*. The inverse of this operation, that associates a quantum operator to a classical 
function, is called 'normal ordering'. In this case one first writes the classical function in 
terms of z and z*, with all the z* 's to the left of the z's, and then replace z by a and z* by 

at. 

The second possibility, called the P representation of A, is obtained by a similar procedure, 
but this time the monomials of A are written in the opposite order, such that they look like 
c nm a n a) m . Once the operator has been put in this form one replaces again a by z and at by 
z* to obtain Ap(q,p). The inverse of this operation is called 'anti-normal ordering'. Notice 
that the differences between the two representations come from the commutator of q and p, 
which is proportional to %. Therefore, these differences go to zero as h goes to zero. 

There is, finally, a third representation which is the most symmetric of all, and therefore 
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the most natural. It is given by the Wigner transformation 



s 



A w (q,p) = j dse*»{q-- A q + -) . (4) 

Aw(q,p) is called the Weyl representation of A j^lO]. Its inverse transformation consists in 
writing the classical function in terms of z and z* considering all possible orderings for each 
monomial and making a symmetric average between all possibilities before replacing z and 
z* by the corresponding operators. As an illustration of these three representations we take 

a 1 9 2 1 2 4 

(m — h — 1) for which we obtain 

Hq = \(p 2 + x 2 ) + x 4 + \(b 2 + b~ 2 ) + 3b 2 x 2 + 36 4 /4 
Hp = \{p 2 + x 2 ) + x 4 - \{b 2 + b- 2 ) - 3b 2 x 2 + 36 4 /4 
H w = l -(p 2 + x 2 ) + x 4 

where b is the width of the coherent state. Notice the term proportional to x 2 that appears 
with opposite signs in Hq and Hp, really modifying the classical dynamics with respect to 
H w . 



B. Basic Path Integrals and their Semiclassical Approximations 

The calculation of the semiclassical propagator in the coherent state representation start- 
ing from path integrals was discussed in detail in 2\. In this section we summarize these 
previous results emphasizing the non-uniqueness of the semiclassical limit as a consequence 
of the overcompleteness of the coherent state representation. The reader is referred to Q] 
for the details. 

In order to write a path integral for K(z" , T; z', 0), the time interval has to be divided into 
a large number of slices and, for each slice, an infinitesimal propagator has to be calculated. 
As pointed out by Klauder and Skagerstam [l|, 0| , there are at least two different ways to do 
that. Each of these gives rise to a different representation of the path integral. Although they 
correspond to identical quantum mechanical quantities, their semiclassical approximations 
are different. We review the construction of these two representations below. 
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The first form of path integral is constructed by breaking the time interval T into N 
parts of size r and inserting the unit operator 

1 = f W^W (8) 



2m 

everywhere between adjacent propagation steps. We denote the real and imaginary parts 
of z by x and y, respectively. In all integrations, dz dz* /2m' means dxdy/m After the 
insertions, the propagator becomes a 2(N — l)-fold integral over the whole phase space 

K(z",t;z',0) = {U^r} nfa^-**^)} ( 6 ) 

J j=l j=0 

with zjy = z" and Zq = z'. Using the coherent state overlap formula 



and expanding e~ tHT l h ~ 1 — iHr/% we write 



(7) 



(z j+1 \e x H M T \ Zj ) = exp j ^(z* +1 - z*) Zj - -z] +1 (z 3+1 - z 3 ) - "jrH j+1J j (8) 
where 

_ (zj+ilHjtj)^) _ 
Hj+ij = — 7- — r-T — = n{z j+1 , zj, tj) (V) 

[Z j+ i\Zj) 

and (1 — iTCj+ijr/fi) has been approximated again by e~ %Hi+1 > : > T ' h . With these manipulations 
the first form of the propagator, that we shall call Kq, becomes 

N-l j_ j_* (N-l r , n 



dzj dz* 
2vr2 

3=1 3=0 



2(^+1 2^+ 1 ^ 2: - 7+1 ^ h^ j+1 ' j 

(10) 

When the limit iV — > 00 and r — > is taken, the above summations turn into integrals. 
Also, 7~Cj+i,j turns into the smooth Hamiltonian function TC(z,z*) = Using the 

properties a|z) = z\z) and (^|a^ = (z\z* , we see that 7i can be easily calculated if H is 



written in terms of creation and annihilation operators with all a) 's to the le 



t of the a 's. 



Therefore, TC is exactly Hq(z,z*), the Q symbol of the Hamiltonian operator 

The second form of path integral starts from the "diagonal representation" of the hamil- 
tonian operator, namely 

H= j \z)h(z\z)^f{z\. (11) 



Assuming that H is either a polynomial in p and q or a converging sequence of such poly- 
nomials, this diagonal representation always exists. The calculation of h is not as direct 

Q 

as that of 7i, but it can be shown |7| that h(z*, z) is exactly Hp, the P symbol of H . To 
facilitate the comparison between this form of path integral, that we call Kp, and Kq, it is 
convenient to break the time interval T into iV — 1 intervals, rather than N. We write 

JV-l 

K P (z",T-z',0) = (z"\l[e-^\z') (12) 

3=1 

and, following Klauder and Skagerstam, we write the infinitesimal propagators as 



/ / /; \ dzi dz* r i-rui * xdzjdz* 

' 1 - / \zj) 1 - -h(z*,z,) -L^-( Zj \ ~ \ \ Zj )e-^^^L {zj \ . (13) 



The complete propagator Kp becomes 

N-l 



rf/dzjdz? r ir } 

K P (z N ,T;z o ,0) = / 11 27ii (zj+i\zj) exp < -—h{z*,Zj) > 

(z* j+l - Zj)zj - -z* +1 (z j+1 - zj) - -^H z p Z j) 



A' -J i _ j_* (N-l .. . 

1 , . 1 ^ , %T . 



j=\ V j=0 



(14) 



Notice that while the two arguments of Hq in Kq belong to two adjacent times in the 
mesh, the two arguments of Hp in Kp belong to the same time. Although both forms 
should give identical results when computed exactly, the differences between the two are 



assical 

sec 



2] 



important for the stationary exponent approximation, resulting in different semic 
propagators. The semiclassical evaluation of Kq and Kp were presented in detail in 
also [9|, llC| |ll|). Here we only list the results: 

K Q (z",t;z',0) = J2]fl^ ^{^S Qu + I Qu )-\{\zr + \A 2 )} , (15) 

K P (z", t; z>, 0) = £ exp { \{S Pv - I Pv ) - \ (|,'f + \zf ) } , (16) 

where 

S lv = S iv (v", u, t) = J dt' [f{iiv - vu) - H % {u, v, t')] - f{u"v" + u'v') (17) 
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is the action and 

i r T 9 2 Hi , , , 

h = 2 1 a^ dt < 18 > 

is a correction to the action. The index i assume the values Q and P and sum over v rep- 
resents the sum over all 'contributing' (complex) classical trajectories satisfying Hamilton's 
equations 

dHi 

ihu = H — — 
ov 

ihv = — 19 
ou 



with boundary conditions 

u(0) = ^ = «' , i>(t) = = v" . (20) 

The factors are an important part of the formulas and they are absolutely necessary to 
recover the exact propagator for quadratic Hamiltonians. If one neglects it, even the Har- 
monic oscillator comes out wrong. For a discussion about contributing and non-contributing 
trajectories, see refs. QQ- 

Finally we remember that the Weyl Hamiltonian can be obtained from H by completely 
symmetrizing the creation and annihilation operators. It turns out to be an exact average 
between Hq and Hp if H contains up to cubic monomials in a and a\ but only an approx- 
imate average for other cases. The semiclassical formula with Hq comes with a correction 
+Iq and that with Hp comes with a correction of —Ip. This suggests a third type of semi- 
classical approximation for the propagator, where one uses the Weyl Hamiltonian and no 
correction term, since the average of +I\ and —To should be approximately zero. This is the 
Weyl approximation, which was conjectured in |2j: 



V v J 

with Sw given by Eq. ()17j1 with Hi replaced by Hw- 

Of the three semiclassical approximations presented, the Weyl approximation seems to 
be the most natural, since it involves the classical hamiltonian directly and no corrections to 
the action. However, this formula does not follow from the two most natural forms of path 
integral proposed by Klauder and used in this section. In the next section we propose a third 
form of path integral which is constructed directly in terms of Hw and whose semiclassical 
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limit is indeed the formula above. For a direct comparison between these semiclassical 
formulas for short propagation times see 



III. COHERENT STATE PATH INTEGRALS WITH THE WEYL SYMBOL 

The new form of path integral we describe in this section is based on an expansion of 
the Hamiltonian in a continuous basis of reflection operators R x whose coefficients H(x) 
are exactly the Weyl symbol of H. We first review the algebra of reflection and translation 
operators in quantum mechanics |4], following closely the presentation in ref. 0|. We then 
use these results to construct the path integral. 

A. Translation and Reflection Operators 

Consider the family of translation operators 

where £ = (q,p) is a point in phase space. It can be shown that the form a complete 
basis, in the sense that any operator A can be expressed as 

A = jm^- (23) 

The Fourier transform of the operators form a complementary family of reflection oper- 
ators R x which also form a basis: 

4 = ^-/d£e^4 (24) 

where x = (Q, P) and iA( = Pq — Qp. In terms of these operators we may write 

A = J ^A(x)R x = ^-Lj J d^dxA(x) et^ 4 (25) 

When this expression is inverted to write that A(x) in terms of A, we find precisely the Weyl 
representation, as given by Eq.fjlJ). This is shown in Appendix A. 

It is convenient to write some of these expressions in terms of a, a*, z and z* instead of 
q, p, q and p. We find that 

T € = e ( * aWa) = e" at e~ z " a e" |2|2/2 , (26) 



which we recognize as the displacement operator [15I Q, 0] frequently used in quantum 
optics. Also 

(zklf^) = e *4-*-**-i-\*\ 2 /* (zklzk-!) (27) 

and 



(z k \A\z k ^) = J dxA{x)(z k \z k ^) j ^(^0^-^-1*1 I 2 . (28) 

Since the integral over £ = (q,p) is quadratic, it can be done immediately. Defining 

(the index k is added for later convenience) we find 

(z k \A\z k ^) = 2 [ dWkdw * k A(w k , wt) e -2|« fc | a +^+a* fc -i»]E-l* fc l a /a-l* fc -il a /2-*):* fc -i ( 30 ) 
./ 2m 



where 

du> fc du>£ dQdP 



(31) 



2ni 2nh 

As the notation suggests, this expression will be our starting point to construct the path 
integral. When A is replaced by the infinitesimal propagator Q- %Hr l h ~ 1 — iHr/% and a 
sequence of these matrix elements are multiplied together, we will find that the all the z k s 
and z^'s appear only in quadratic forms and can be integrated over. The resulting path 
integral will be written in the new variables w. 



B. The Path Integral 



We start from 



K(z",t; z',0) 



2ni 



}n{^~^ )T i^> 



(32) 



i=i 3=1 

where z^ = z" , z$ = z' , r is the time step, Nt = T and we take N to be even for convenience. 
The infinitesimal propagators can be calculated with Eq. (J30|) by simply replacing A(x k ) by 
e -iH(x k )r/h w ]j ere H{x) is the Weyl symbol of H calculated at (Q k , P k ). We obtain 



k [f, t, o) = 2" }{nu i^£} /{nt 1 ^} 



x 



(33) 



exp {EL ~^fcT - 2K| 2 + 2z* Wfc + 2* fc _ 1 < -tf- ^ 
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where = H(wk,wl). The integrals over the z/s and the z*'s can be performed exactly. 
When this is done we find 



K(z»,t;z>,0) = S{X\U^}^-^-^ 



(34) 




where 



<Pn = TtLi [-iTH h /h - 2K| 2 + 2z"*w N+1 - k (-l) k+1 + 2z'w* k (-l) k+1 ] 



(35) 



+ z'z' 



In the second line of (|34|) we have written the dependence of the propagator on z' and z' 
explicitly and defined 



C. Alternative form and the limit of continuum 

Eqs.(jHIJ) an d (|33|) correspond to the coherent state path integral in the Weyl representa- 
tion. It is very different from the previous forms presented in section II in two respects: the 
measure lacks a factor 2 in the denominator and, more importantly, the exponent does not 
resemble an action at all. Although these expressions appear to be the most practical for 
actual calculations, we can manipulate the terms in (p^ to make it look more familiar and 
similar to an action function. However, it is only when we take the limit of the continuum 
that we really recognize the action as part of the exponent. We shall do these manipulations 
now, but we insist that Eqs. (j3^jl and (|33|) are the direct analogs of Eqs. (fTU|) and (fT4"|) for the 
Q and P representations respectively. Although unusual, and perhaps more complicated, we 
shall see that, in the semiclassical limit, the Weyl form becomes the simplest of them all. 

We show in Appendix B that the quadratic terms in (p^ can be written as 




Ef=i H^A-2| Wfc | 2 ]+4Ef= 1 1 E-=i^ + i^ + i- J (-i) J 
= EL«*+i-*(-i)* +1 - 



(36) 



: 2 Ef=i - O - w i+i( w k+i - w k )] 



(37) 
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where the sums on the right go in steps of two. The terms proportional to z' and z"* can 
also be re-written as 



(38) 



(39) 



E fc =i w w+ i-fc(-i) fc+1 = E fc =i ! 3(« ; fc+i - 

When these terms are replaced in the exponent we get 

<t>N = 2 J2k=l,3 [ w k(K+l ~ w l) ~ W* k +A w k+1 - W k )] - f £f=l 
- 4 Lfc=l,3K+l - 2^/=fe+l,it+3(^+2 

This is the alternative discrete version of 4>n- Although not much enlightening than the 
original form, Eq.(|35|). the first line shows a closer resemblance to the usual action function. 
More importantly, this expression is ready for the continuum limit. Taking iV — > oo, r — > 
with Nt = T we obtain 

rT tt i , , rT / . * * • \ i , rT • / ,\ f-T 



(40) 



= -f / Hdt + J*q (to* - w*w) dt - J Q «;(*) / t w*(t')dt'dt 
-z' J T w*dt + z"* J T wdt + z'z"* . 

Notice that the factors of 2 and 4 compensate for the sums in steps of two. 
The integrals in the last term on the first line can be rewritten as 

w(t)[w*(T)-w*(t)]dt = w*(T)[w{T)-w{0))- / w(t)w*{t)dt. (41) 
'o Jo 

The last term above cancels one of the terms in Eq. (j4()j) . After performing the integrals on 

the second line of Eq. (j40j) . making some simple rearrangements and an integration by parts, 

we can write the exponent in the form 



<P= j;S+(z' -w(0)) w*(0) + 



z"*-w*{T) 



+ {z"*-w*{T)) w(T) + 



z'-w(0) 
2 



(42) 



where S is the (complex) action 



121 



s 



t r 



ih 



(w*w — ww*) — H 



dt - j(w*(T)w(T) + w*(0)w(0)). 



(43) 



Notice that the action in the coherent state representation is not just the integral corre- 
sponding to pq — H, but it includes important boundary terms. Besides, the exponent <fi of 
the path integral is not just the action and also includes further boundary terms. We shall 
see, however, that the extra terms in Eq. (|42jl vanish in the semiclassical limit. 
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IV. SEMICLASSICAL LIMIT 



The semiclassical limit of the propagator is obtained by performing the integrals over Wk 
and w* k with the stationary phase approximation. Because the exponent <pN is not a phase, 
but a complex quantity, we use the terminology 'stationary exponent approximation'. 



A. The Stationary Exponent Condition 



Using Eq. (|H5jl for N even and I 7^ 1 even we obtain 
d(j) N ir dH[ 



and 



it dH u 



- 2wi - 2z + 4[l0;_i - W1-2 H W2 + W1] =0 



/_1 - 2w t+1 + 2z + A[wi - wi-x H h w 2 — wx] = 0. 

Adding these two equations we obtain simply 



h dwi 



1 1 



For I = 1 we get 



T 



(44) 



- 2wi + 2z = 0. 



dw\ h dw\ 

For the derivatives with respect to wi we proceed in the same way. For I odd we get 
iTdHi 



(45) 



>N 



and 



dwi 



h dwi 



- 2w* - 2z"* + 4K +1 - w* l+2 + ■■■- w* N _, + w* N ) = 



Adding the two equations we obtain 



- 2w* l+l + 2z"* + 4K +2 - + • • • + w* N _ x - w* N ] = 0. 



Finally for I = N we get 



z 1 
ft~2 



dwi dwi+i 

IT dHi 



2u^ + 2z" 



0. 



(46) 



(47) 



Own ft 9wn 

Taking the continuum limit and using the u and v variables in the place of w and w*, 
Eqs.flEj), flUJ), (|SJ) and flSJ) become 

dHyy . dHy\r 



ihu = + 



dv 



ihv 



du 



(48) 
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with boundary conditions 

u(0) = z , v(T) = z"* . (49) 

The average of the derivatives at consecutive time steps that appears in the left side of 
equations (|44|) and resemble the stationary conditions obtained in |3j. In that case, 
however, one of the derivatives involved Hp and the other Hq. 



B. Expansion Around the Stationary Trajectory 

Let w® and w* k ° represent the stationary trajectory and + and w* k Q + £1 a nearby 
path. Expanding the exponent up to second order around the stationary trajectory we get 

<P N = <P N + 5 2 <i )N + O(3) (50) 

(the first order term is zero) with 

pfa = EL {-U A ^l + 2C ^a + s k a 2 ] - 2^a} + 4 e^i 1 c k+1 e U €jH-w(-i) i+1 



= -lx T K N x 



2' 



(51) 



where X T = (£ N , £* N , £ N _ h . . . , £i, and 



A * ~ eSf ' u * ~ a^f ' ° fc _ dw k dwt ^ oz ) 



are calculated at the stationary trajectory. 

When the limit of the continuum is taken, the boundary conditions Eq. (J49|) kill the extra 
terms in the exponent 0, Eq. (j42j) . which becomes simply the action of the complex trajectory. 
Therefore the semiclassical propagator becomes 



K w (z',z",T) = e i^(l*T+l*"l 2 ) lim 2 (53) 



As usual, the calculation of the determinant of the quadratic form is the most lengthy step 
of the semiclassical calculation. In this case the calculation is particularly tricky, because of 
the double sum in the last term of the first line of equation (jB*T|) . To avoid losing the focus 
with this lengthy algebra here we do the calculation in the Appendix C. The final result is 
indeed the conjectured formula, Eq. (j21j) that we repeat here: 
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Of course, if there is more than one stationary trajectory, one should sum over all the 
contributing ones. 



V. A COMPARISON BETWEEN THE THREE FORMS OF PATH INTEGRAL 

In principle, all discrete forms of path integrals given by Eqs. fjlOj) . (j!4j) and (|H4jl are 
quantum mechanically equivalent. For fixed N, however, they are not identical and in the 
limit N — > oo there are well known convergence problems, making the comparison difficult. 
In order to illustrate the differences between the three forms we shall study the discrete 
propagators for the simple harmonic oscillator. The Hamiltonian operator is 

h 2 d 2 muj 2 x 2 . ( ± 1 

H = —r H = fvuj [aJa + - 

2m dx 2 2 V 2 

and, choosing the coherent state width as b = y^h/moo, the classical symbols, in the u and 
v variables, are 

Hq = Too ^uv + Hp = Tiuj (^uv — -J H w = Tiujuv. 

Using H w in the stationary conditions (f4*4" j) - (}4"Tj) we obtain the stationary path 

*k-l fy*N-k 

Wk = ~^F~ z Wk = z 

where 

a = 1 + iruo/2. 

The calculation of the phase cjP N at the stationary trajectory is lengthy but involves only 
simple geometric sums. Several simplifications occur when all the terms in <fi° N are added 
together and the result is 



^) z'z"*-l(\zf + \zr). 



The determinant of the quadratic form is calculated in Appendix C and results in (see 
Eqs.flCI]) and (fHHj) ) 

det A N = 2 2N i 2N a 2N . (55) 
Putting everything together we obtain 

ts- I H I rp\ /, , • /o^AT (±=i^) N z'z''*-\z'\ 2 /2-\z"\ 2 /2 

Kw{z , z ,T) = (1 + ituj/2) e \i+^/2j (56) 
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which clearly converges to the exact propagator as r — > 0. Doing similar calculations for 
the Q and P propagators we find 

Kq(z", z', T) = e -i-T,2Hl-ir^ Z '^-\AV2-WT/2 (57) 

and 

K P (z", z>, T) = (1 + iru)- N e ^ T/2+ (^- |2 ' |2/2Hz " |2/2 (58) 

which also converge to the exact result. Notice that the overall phase —iuT/2 comes out 
exact for Kq even in the discrete form. However, the term multiplying z'z"*, which goes to 
e~ lujT as N — >• oo, converges much faster for K\y then the corresponding terms in Kq or Kp. 
Moreover, for any finite value of N, this term has unit modulus in Kw, while its modulus 
is larger than one for Kq and smaller than one in Kp. Just for the sake of comparison 
let us call fi this coefficient. Taking uoT = 2tt and N = 100 we find \xq ~ 1.22 + O.Oli, 
Hp « 0.82 + 0.007i and \i w « 0.999998 + 0.002i This suggests that the new path integral 
representation should be better than the two KS forms for numerical evaluations. 



VI. CONNECTING THE WIGNER AND THE HUSIMI PROPAGATORS 

In this section we show that the Weyl representation of the evolution operator 

U{q,p,T) = J(q- s/2\e- iflT/h \q + s/2) e ips/h ds (59) 

can be directly related to the path integral representation derived in the section III. This is 
an interesting formal result that was also obtained by Ozorio de Almeida in section 6 of 
starting from the opposite direction, i.e., from the path integral representation of U. The 
result provides an explicit connection between these two famous phase space representations 
of quantum mechanics. As we shall see, the connection is very simple when written in terms 
of path integrals. 

We start by re- writing Eq. fJoTJj) as 



U(q,p,T) = J(q-s/2\z")(z"\U\z')(z'\q + s/2) e^ds^ 



'■V* &z"&z"* 



2m 2ni 



JV[w,w*}e^J dse ips / n J ^f- ^gf- (/3 1 z" } (z'\a) x 
exp 



(60) 



Jjf _ kjp + z > z »* + 2Cz"* - 2C*z' 
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where we used Eqs. ljMj) and (|3l)|) and denned a = q + s/2 and /3 = q — s/2 in the second 
line. The integrals in z' and z" are quadratic and be performed analytically. The integral 
over z' is straightforward and gives 



■ X 



(P\z")exp[-^ + 2Cz"*-^-^^ + ^(z"*-2C*)j. 
It can be seen by inspection that the exponent in the second line above can be written as 

n 1 'W 2 (z"\a + A)e- B (62) 

with A = bV2(C + C*) and B = -A 2 /2b 2 - Aa/b 2 + 2C* 2 + 2y/2aC* /b. When (E2D is 
substituted into ([6ip the integral in z" produces (/3\a + A) = 6 (a — j3 + A) = 8(s + A). 
The delta function takes care of the integral over s and after some simplifications we obtain 
simply 

U{q,p,T) = J V\wMY^ 2Cz% - 2C * z ^ C? (63) 

where z x = (q/b + ipb/h)/y/2. Comparing with Eq.(j2U) shows that the path integral for 
U(q,p,T) is directly related to that for K(z x , z x ,T). Indeed, the path integral for U has 
a single extra term 2\C\ 2 with respect to the K. One might say that this terms promotes 
the 'unsmoothing' of the coherent state propagator. Conversely, the diagonal coherent state 
propagator has the extra term — 2\C\ 2 with respect to the U, smoothing it out. This result 
was also obtained by Ozorio de Almeida in section 6 of [5]. The coefficient C can actually 
be interpreted as the Wigner chord linking the ends of a polygon in phase space whose sides 
are given by the Qk and variables in Wk- We can also calculate explicitly the two terms 
that involve q and p in (|63|) . Using the definition of C in Eq. (J36|) we find that 

N 

Ti 



2Cz* x - 2C*z x = ^{QkP ~ P k q) (64) 



k=l 

which is the sum of the symplectic areas between X k = (Qk,Pk) an d x = (q,p) and is 
independent of the width b. 



17 



APPENDIX A: EXPANSION IN REFLECTION AND TRANSLATION OPER- 
ATORS 



This appendix follows closely the demonstration in |5J. A comparison between Eqs. 
and shows that 

A(0 = ^ J dxA(x)e* x/ * (Al) 
and, inverting the Fourier transform, 

M*) = A / d^(Oe-i^. (A2) 



2tt/i „ 

Using Eq. ()23)1 again in the coordinate representation we obtain 

= J &£A{q, p)S(q + -q.-q) e M«-+* ) (A3) 

This Fourier transform can be inverted as follows: we define q' = q + — g_, Q = (q + + g_)/2, 
multiply both sides by e~ tp 'Q/ h and integrate over Q. The integral over Q on the right-hand- 
side yields a delta function on p — p' and we obtain 

A(0 = JdQ(Q + q/2\A\Q - g/2>e-*rf . (A4) 

where (q',p') has been changed back to Finally we use Eq. (|A2j) to get 

= ^ / dgdpdQ e iKQ-0)-iP 9 (q + q/2 \A\Q - q/2) 
= fdqe-* p * (Q + q/2\A\Q - q/2) 

which is the same as Eq.ffl. 
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APPENDIX B: PROOF OF EQ. (f57j) 



First we re-write, for N even, 

4w 2 W 1 + 

4w 3 [w 2 - Wi] + 

Awl[w 3 - (w 2 -w 1 )} + 

4wZ [{w A - w 3 ) + (w 2 - i«i)] + 

4w*[w 5 - (to 4 - iu 3 ) - (iu 2 - u>i)]+ 

iw* N [w N -i - (W N _ 2 ~ WN-3) (W 2 - Wl)) = 

A[w 2 Wi + W* 4 W 3 + WgW 5 H h W* N W N -i] - 

A(w 2 - Wl ) [(wl - w*) + (w* 6 - w* 5 ) + • • • + (w* N - w* N _ ± )\ - 
4(w 4 - w 3 ) [« - wl) + {w* 8 -w* 7 ) + ■■■ + (w* N - w* N _i)] - 



4 Ef=i w t+i w k - 4 Ef=i 3 3 K+i - w k) EjI fc +i,fc+3«+2 - <fi) 

The second term is already in the form needed for Eq. (j37j) . The first term is now modified 
as follows: half of it remains unchanged and, in the second half, we add and subtract terms 
as in 

w* k+1 w k = w* k+2 w k+1 - [w k+1 (w* k+2 - w* k+1 ) + w* k+1 (w k+1 - w* k )} (B2) 
for k = 1, 3, . . . , iV — 3 only. We obtain 

4 L fc =i,3 w k+iWk = 2 ^ fe=13 w k+1 w k + 2 ^ fc=1>3 w k+2 w k+1 - ^ 

2 Ef="l 3 3K+l«.+2 - <+l) + <+l(™k+l ~ w* k )] + 2vf N w N . x 

We finally add — 2J2k=i w t w k- The part of this sum containing odd k's goes together with 
the first sum above. The even k's up to N — 2 goes with the second sum. We get 

4 EEs w t+i w k - 2 Ef=i w >k = 

2 Ef=i M™t+i -w* k )+2 Ef = t 3 w k+1 (w* k+2 - w* k+1 )- (B4) 
2 Ef=i! 3 K+i«+2 - + - «>*)] - 2w* N (w N - W N ^i) 
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The second term in the second line cancels against the first term of the third line. After 
incorporating the last term into the sum we get 

4 EEs w t+i w k - 2 Ef=i w* k w k = 2 EEsHK+i - K) - w* k+x {w k+1 - w k )\ (B5) 

APPENDIX C: CALCULATION OF THE DETERMINANT 

The quadratic form in Eq. (j53|) is defined by the matrix 

f irA N /h irC N /h + 2 ... X 

irC N /h + 2 irB N /h -4 4 -4 

-4 irA N ^i/h irC N -x/n + 2 

irC N ^/h + 2 irB^h -4 4 

4 -4 




-4 
... -4 irAx/h ird/h + 2 
\ ... ird/h + 2 irBx/h J 

whose determinant, det An, we seek. To simplify the notation we will drop the det symbol 
in this appendix and use simply A^v for det A at. It is useful to factor 2i out of each element 
and call the new determinant A^v- Of course 

A N = 2 2N i 2N A N . (CI) 

This cancels both the 2 N and the sign (— 1)^ in Eq. (j53J) . leaving only A at- Next we do the 
following sequence of operations that do not change the value of the determinant: 
column 2 — > column 2 + column 4 
column 4 — > column 4 + column 6 



column N-2 — > column N-2 + column N-4 
line 2 — > line 2 + line 4 
line 4 — > line 4 + line 6 
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line N-2 -> line N-2 + line N-4. 
This put the matrix in block tri-diagonal form: 

/ rJ t~^ n n n \ 

tC n _ ■ t(B jy+Bjy„i) tCn-i I • t_Bjv-i 

2ft 2ft 2h 2ft 

pi tCjv-i i • t Ajy_i rCjy_i 

U 2% 1 2% 2h 

n T -Bjv-1 tCn-i • " 



2h 2h 2h 2h 
































tB]V-2 

2ft 






^ + l 







r Cjy, 

2h 

tB n _ 2 



2% 



TO-l I „• T±tl 

2ft ~ 6 2ft 

tCi I • tAi rCi _ • 

2h 27i 2^ 6 

\ rBi rCi • rgl / 

\ ■ ■ ■ 2% 2h L 2h J 

We can now compute the determinant using Laplace's method. Let Y ^ be the determinant 
obtained from the matrix above by removing the first line and the first column. The two 
determinants An and Tn satisfy the following recursion relation: 

T N = ligsd^ A^.x- (l^i + j) 2 r jV _ 1 + (C2) 



"" 2 Cjv-l 1 \ tB n _ 1 a i tB n _ 1 



1 + ^(^jV-l - 



AjV-2 



Keeping only terms of first order in r and taking the limit r — > we find 

^"f"- 1 = ^r* + i^A^.x + G(r 2 ) 

r^-r^ = (b^+b^o Ajy _ t _ + ^A^ + 0(r 2 ) 

or 

A =ir + ^A 

with initial conditions A(0) = 1 and T(0) = 0. 

Notice that in the case of the harmonic oscillator Hk = TiuJWkW* k and, therefore, Ak = 
13^ = and = hu. In this case Eqs. flC2| ) can be solved exactly, without the need to take 
the continuum limit. We find simply 



(C3) 



(C4) 



IN 



21 



2N 

VjJT 



which can be iterated to give 

A N =(l + ^j • (C6) 

To solve Eqs. ()C4|) in the general case we need a last change of variables Q = 2iA. In the 
new variable we get 

% h / C7 x 
h h 

with fi(0) = 2z and T(0) = 0. Identifying T with u and with t>, we recognize these 
equations immediately as the equations of motion (|48)l linearized around the stationary 
trajectory The solution we seek, A(T) = Q(T)/2i can be obtained with the help of the 
relations 

where we use a single prime for quantities calculated at t — and a double prime when 
t = T. A variation in the second of these equations leads to 

Using 5u' = T(0) = 0, 5v" = fi(T) and 5v' = fi(0) = 2i we get 
and 
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